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0 Method and apparatus for the acoustic investigation of a casing cemented in a borehole. 



@ Method for determining at least one characteris- 
tic of a casing cemented in a borehole, such as 
cement bond and casing thickness, from a reflected 
acoustic signal S(t) obtained by directing ^n acoustic 
pulse at a substantially normal incidence towards a 
radial sector of the wall of the casing, said pulse 
stimulating thickness resonance within the walls of 
the casing. 

The method according to the invention com- 
prises the steps of : defining a first time window 
^corresponding to a first portion of signal S(t) includ- 
^ing the Initial reflection from the casing and subse- 
g^quent acoustic returns due to resonance, defining a 
second time window corresponding to a second por- 
^tion of signal S(t) only including the initial reflection 
l^from the casing, and determining said casing char- 
ge acteristic from information related to resonance con- 
POtalned in said first time window while normalizing 
Q^said information by information contained within said 
second time window. 
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METHOD AND APPARATUS FOR THE ACOUSTIC INVESTIGATION OF A CASING CEMENTED IN A 

BOREHOLE 



The invention relates to the acoustic Inspection 
of a casing cemented in a borehole which traverses 
underground strata. It relates more particularly to a 
method for processing a measured acoustic signal 
for forming signals representative of characteristics 
relating to the casing and to the sunrounding ma- 
terials, while simultaneously compensating for vari- 
ations in the acoustic pulse generated and the 
acoustic attenuation by the mud filling the 
borehole. 

In general, once a well has reached the desired 
depth, the borehole is cased, with cement being 
injected into the annular space between the casing 
and the wall of the borehole in order to prevent 
hydraulic communication between the various geo- 
logical strata. In order to determine whether such 
unwanted communication nevertheless ejcists. mea- 
surements may be performed downhole by means 
of a logging tool in order to determine the quality 
of the bond between the cement and the casing. 
Another information of interest which can be de- 
rived from these measurements is the thickness of 
the casing- 
It has long been the practice to use acoustic 
waves for performing such measurements. Such 
techniques have relied on measurements which are 
averages in the circumferential and/or longitudinal 
direction of the casing, and consequently cannot 
identify localized phenomena such as longitudinal 
hydraulic communication paths. A complete pre- 
sentation of prior art techniques can be found in 
U.S. Patent 4.255.798 to Havira, columns 2-7. 

Of the techniques that have sought to improve 
the vertical and radial resolution in such inspection, 
the technique described in the above-mentioned 
Havira patent turns out to have been an extremely 
important breakthrough. This technique consists in 
emitting an acoustic pulse over a radial sector of 
the casing, with the pulse being constituted by 
acoustic waves at frequencies selected to cause 
resonance to appear across the thickness of the 
casing between its outside wall and. its inside wall; 
in determining the energy present in a reverbera- 
tion segment of the reflected signal; and in char- 
acterizing from said energy the quality of the bond 
of the cement behind said radial sector of the 
casing. The reverberation segment under consider- 
ation is selected so as to be substantially repre- 
sentative of acoustic reverberation between the 
walls of . the casing. Rapid damping of the reso- 
nance, in other words low energy, means that ce- 
ment is present behind the casing, whereas slow 
damping, i.e. high energy, means that there is no 
cement. 



A logging tool using such a technique is de- 
scribed in a commercial brochure entitled "Cement 
Evaluation Tool" published by Schlumberger in 
June 1983. The sonde of this tool is centered 

5 within the casing. 

The sonde includes eight transducers distrib- 
uted helically at 45* intervals, thereby obtaining 
good coverage around the periphery of the casing. 
Acoustic pulses are fired sequentially. They are 

10 likewise received sequentially and analyzed and 
transmitted to the surface where they are pro- 
cessed. 

In addition, a ninth transducer referred to as 
the "reference" transducer points along the axis of 
IS the casing towards a reflecting wall which is plane 
and disposed at a fixed distance from the reference 
transducer. The reflection signal detected by the 
ninth transducer is used to determine in situ the 
propagation time through the borehole fluid (mud). 
20 i.e. the time interval between emission and recep- 
tion of the acoustic wave, and the wave's propaga- 
tion velocity through the mud is deduced there- 
from. Given the propagation velocity of the wave, it 
is possible to deteonine the apparent radius of the 
25 casing for each of the eight transducers. It is par- 
ticularly advantageous to obtain this radius since it 
makes it possible, in particular, to detect any possi- 
ble deformation of the casing and to monitor the 
centering of the sonde inside the casing in order to 
30 obtain an indication of the validity and the quality of 
the measurements as performed. 

In this technique, a portion is taken from the 
reflected acoustic signal, which portion is repre- 
sentative of acoustic reflections between the wails 
35 of the casing. 

A signal QCi is derived therefrom representa- 
tive of the quality of the bond of the cement with 
the casing on the basis of the energy W2 mea- 
sured in a reverberation segment of the reflected 
40 signal S. 

The changes in amplitude of such a signal S 
as a function of time and as picked up by a 
transducer are shown in Figure 2. The effects due 
to the mud are taken into account by normalizing 
45 the measured energy W2 relative to the peak am- 
plitude signal Wl of the reflected acoustic signal. 
Nevertheless, it turns out that this nonmalization 
does not give complete satisfaction since quantita- 
tive interpretation of the measurements shows up 
so problems of divergence in the measurements. 

Tests and experiments performed by the Ap- 
plicant have shown that these differences can be 
attributed to a large extent to the fact that the 
conventional processing is sensitive to the prop- 
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erties of the mud. to the characteristics of the 
transducer and of its drive electronics. 

The Applicant has observed that the normaliza- 
tion of the portion of the signal at energy W2 was 
being performed using a peak amplitude signal W1 
which did not correspond to the same frequencies 
as those present in the portion of the signal cor- 
responding to the energy W2. 

In general, the energy W1 is not at the same 
frequency as the energy W2. even in water. The 
energy W1 is a measure of the energy maximum 
conveyed by a spectrum component which de- 
pends on the properties of the pulse i.e. the char- 
acteristics of the transducer and its drive elctronics. 
and on the attenuation by the mud. whereas the 
energy W2 includes energy only around the reso- 
nance frequency of the casing. 

The Applicant proposes a different concept for 
processing the reflected acoustic signal, whereby 
variations in acoustic pulse characteristics and mud 
attenuation can be compensated for in a satisfac- 
tory fashion and the reliability of the cement bond 
and other information obtained from the measure- 
ments can be improved. 

There is provided according to the invention a 
method for determining at least one characteristic 
of a casing cemented in a borehole, such as ce- 
ment bond and casing thickness, from a reflected 
acoustic signal S(t) obtained by directing an acous- 
tic pulse at a substantially normal incidence to- 
wards a radial sector of the wall of the casing, said 
pulse stimulating thickness resonance within the 
walls of the casing, comprising the steps of : defin- 
ing a first time window corresponding to a first 
portion of signal S(t) including the initial reflection 
from the casing and subsequent acoustic returns 
due to resonance, defining a second time window 
corresponding to a second portion of signal S(t) 
only including the initial reflection from the casing, 
and determining said casing characteristic from 
Information related to resonance contained in said 
first time window while normalizing said information 
by information contained within said second time 
window. 

Other characteristics and advantages of the 
invention appear from reading the following de- 
scription made with reference to the accompanying 
drawings, in which: 

Figure 1A is a diagram of an apparatus for 
cement evaluation, the apparatus being shown in 
use in a borehole; 

Figure IB shows in more detail the sonde of 
the apparatus of Figure 1A; 

Figure 2 shows a reflected acoustic signal S 
as detected by the transducer T. with signal por- 
tions W1. W2 as defined in the prior art method: 

Figure 3 shows a reflected acoustic signal S 
as detected by the transducer T, with signal por- 



tions CW and PW defined in accordance with the 
present invention; 

Figures 4 and 5 show, respectively, the nor- 
malization window C and the processing window P 
5 as defined in the method of the invention; 

Figure 6 shows the frequency spectrum of 
the signal amplitude, with solid lines for the signal 
windowed with the processing window and with 
dashed lines for the signal windowed with the nor- 
10 malization window; f 

Figure 7 shows a portion of the group delay 
spectrum of the signal after normalization, high- 
lighting a dip indicative of a resonance; 

Figures 8 and 9 show in solid lines the 
;s amplitude spectrum and in dashed lines the group 
delay spectrum of the reflected signal as obtained 
from synthetic data for, respectively, casings of two 
different thicknesses; 

Figure 10 is a diagram of an electromagnetic 
20 model of the casing, the cement, and the mud. for 
use in a method for computing the cement imped- 
ance and thickness of the casing. 

With reference to Figure 1 . an acoustic logging 
apparatus for studying the quality of the cement 
25 bond comprises a downhole device 10 suspended 
in a borehole 11 at the end of a multiconductor 
cable 12. 

The borehole, filled with a fluid, typically drill- 
ing mud M. passes through strata 13 and is lined 

30 with casing 14. Cement 15 occupies the annular 
space between the casing and the geological strata 
traversed by the borehole. This logging apparatus 
is based on the measurement concept of casing 
resonance as disclosed In the above-mentioned 

35 Havira patent. 

In general, the downhole device 10 is an elon- 
gate body including a sonde 16, top and bottom 
centralizers 17 and 18 for holding the sonde 16 
centered on the axis of the casing 14, an acoustic 

40 compensation section 22. a hydraulic compensa- 
tion section 32, and an electronics section 19. 

•Various embodiments are possible for the 
sonde 16. In a first embodiment (not shown) the 
sonde may comprise eight acoustic transducers 

45 angularly distributed around the sonde and dis- 
posed helically so as to enable eight sectors of 
casing to be inspected, as described in the afore- 
mentioned Havira patent or U.S. patent 4.685.092 
to Dumont. 

50 In a second embodiment, the sonde 16 has a 

single acoustic transducer T. The sonde 16 has a 
drive motor 20 for rotating the assembly carrying 
the transducer T about the axis of the sonde. 

ThuSt in normal operation, the sonde can in- 
55 vestigate the entire periphery of the borehole by 
continuously rotating the transducer. 

A more detailed description of the structure of 
the sonde in this embodiment is given in US ap- 
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plication N* 494.721 filed March 16. 1990. By 
having a system with an offset drive shaft and an 
abutment it is possible in a calibration mode of 
operation to decouple the reflector 24 from the 
drive shaft so as to cause only the transducer T to 
be rotated into a position where It looks at the 
reflector. 

The centralizers 18 and 19 are conventional 
and are provided with means (not shown) for estab- 
lishing electrical and/or hydraulic connections be- 
tween the top portion and the bottom portion of 
each of the centralizers. 

The hydraulic compensation section 32 is of 
conventional type and serves to keep the pressure 
of the fluid in the hydraulic circuits in the downhole 
device 10 at the same pressure as the borehole 
fluid M. i.e. the hydrostatic pressure of the well, 
thereby avoiding excessive differential pressures 
on sensitive portions of the sonde. 

The section 22 containing the acoustic trans- 
ducer has a hollow portion 23 which is open to the 
mud M. The transducer T emits acoustic waves 
through the hollow portion 23 along an axis which 
is substantially perpendicular to the axis of the 
sonde 16, either towards the reflector 24 in the 
calibration mode or towards the casing (normal 
measurement mode). The reflector 24 which is 
substantially perpendicular to the emission axis of 
the transducer T simulates a radial sector of the 
casing 14 under investigation. 

The transducer T is controlled to emit acoustic 
pulses of short duration and having a spectrum 
covering a frequency range from 200 kHz to 700 
kHz, with the mean frequency being about 500 
kHz. 

Since the transducer T emits radially, each 
emitted pulse is reflected by the various interfaces 
it encounters, namely the casing 14 or the reflector 
24. thereby giving rise to an echo signal which is 
detected by the same transducer. 

The detected signals are transmitted to the 
electronics section 19 by appropriate cabling (not 
shown). In the electronics section 19 the signals 
are sampled at a selected rate, put into digital 
form, multiplexed, and finally transmitted to the 
surface via a modem (not shown) connected to 
conductors in the cable 12. 

On the surface, the conductor's of the cable 12 
are connected to an acquisition system 25 includ- 
ing a modem for decoding the transmitted informa- 
tion, and also including a demultiplexer. The sig- 
nals are then stored in a memory 26. The memory 
is connected to a computer 27 suitable for perform- 
ing a sequence of processing stages at each 
depth. These stages are explained in greater detail 
below. 

Signals representative of the depth of the de- 
vice 10 are generated by a displacement detector 



comprising a pulse generator 28 controlled by a 
wheel 29 bearing tangentially against the cable 12. 
The pulse generator 28 delivers a pulse each time 
the cable 12 has moved through a unit length. 

5 The processing performed by the computer 27 

provides a first set of output signals which are 
delivered to a magnetic recorder 30 which stores 
them, for each depth, in digital form on a magnetic 
tape. The magnetic recording comprises, in par- 

w ticular, the signals ^ delivered by the transducer T 
and it consequently enables the data to be pro- 
cessed subsequently. Also for each depth, the 
computer 27 provides a second set of output sig- 
nals to an optical recorder 31 which provides a 

15 graphical recording on a film. e.g. of the type 
described in French patent 2 491 123. 

Figure 2 shows the variations as a function of 
time in the amplitude of a signal S as sensed by 
the transducer T. and illustrates the method of 

20 windowing signal S disclosed in the Havira patent 
A signal representative of the quality of the 
bond between the cement and the casing is pro- 
duced from the energy W2 measured in a portion 
of the reflected signal S which essentially cor- 

25 respond to reverberations between the walls of the 
casing as excited by the pulses emitted by the 
transducer. 

In the window ta - ta. the energy W2 may also 
include energy due to multiple reflections on the 

30 geological strata. These reflections can be ob- 
served from instant ti . 

The energy W2 is normalized by the maximum 
peak energy W1 in the sensed signal S. This 
normalization serves to eliminate effects due to the 

35 mud. but it is not fully satisfactory as explained 
hereinabove. 

Figure 3 shows the same variations in am- 
plitude of the signal S as a function of time and as 
sensed by the transducer T. together with two 

40 windows defined as explained below, which both 
include a time reference or time origin to located in 
a portion of the signal which corresponds to the 
initial reflection from the casing wall. 

45 

Determining a time origin 

According to the invention, a stable time origin 
is determined so as to escape both from noise and 
50 from variations in transducer response over dif- 
ferent measurements. 

To do this, several methods can be contem- 
plated. 

In a first possibility, an analytic signal Sa(t) is 
55 determined as follows: 
Sa(t) = S(t) - iH(S(t)) 

where S(t) is the measured signal and H is the 
Hilbert transform of the measured signal, with i 
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being the square root of -1 . The Hilbert transform is 
obtained by any known algorithm. 

The time origin to is defined by the instant at 
which the analytic signal has a peak. 

Another method of obtaining the analytic signal 
Sa(t) consists in performing a Fourier transform by 
any fast Fourier transform algorithm, then in remov- 
ing the negative frequencies by filtering, and sub- 
sequently in performing a further fast Fourier trans- 
form on the signal obtained after filtering. 

The instant at which a peak occurs in the 
analytic signal is taken as being the origin to. 

An even faster method consists in rectifying 
the signal received by the transducer, in extracting 
the amplitude maximum frorn the rectified signal, 
together with two extreme values on either side of 
said maximum so as to obtain three points defined 
by (time, amplitude) coordinates. 

A second order polynomial passing through 
these three points is then determined in conven- 
tional manner. 

The maximum of the resulting polynomial is 
determined and the time coordinate to of the maxi- 
mum is taken as being the time origin and serves 
as a reference for determining the position and 
duration of the windows. 



Windowing the signal S(t) 

A first portion PW of the signal S(t) is defined 
by a wide time window (-ti. ta) designated as 
processing window hereinafter. The signal portion 
PW obtained in this way includes the initial reflec- 
tion of the acoustic signal from .the inner wall of the 
casing and acoustic returns due to reverberations 
between the walls of said casing caused by its 
resonance in response to the acoustic pulse. 

A second portion CW of signal S(t) used for 
nonmalization purposes is defined by a narrow time 
window (-ti. -f-ti). also called the normalization 
window, which window is centered at time origin to 
and can be defined as the early portion of the wide 
window. The signal portion CW thus defined mainly 
comprises the initial reflection from the wall of the 
casing. 

Figures 4 and 5 respectively show an example 
of the windowing signal C(t) corresponding to a 
narrow window as used for normalization, and an 
example of the windowing signal P(t) correspond- 
ing to a wide window as used for performing the 
processing described below. 

Windowed signals corresponding to the por- 
tions CW and PW may be obtained by multiplying 
the signal S(t) by the signals shown in Figures 4 
and 5. namely the narrow window signal C(t) and 
the wide window signal P(t). 

The time intervals spanned by the windows are 



defined as multiples of the nominal resonance pe- 
riod of the casing, so as to automatically adapt the 
windows to the type of casing under investigation. 
The time intervals -ti ,tb and to.ti should lie in 

5 the range between 1.5 and 3 times the nominal 
resonance period of the casing, and are preferably 
equal to 2.5 tlm^is this nominal period. The nor- 
malization window -ti,ti thus lasts for a period lying 
between 3 and 6 times the nominal resonance 

;o period of the caslngi preferably a period equal to 5 
times this nominal period. 

A suitable range for the time interval to.ta of the 
processing window is from 6 to 12 times, prefer- 
ably 6 to 8 times the nominal resonance period of 

rs the casing. Thus, the entire time interval -ti.ta 
spanned by the processing window is suitably be- 
tween 7.5 times and 15 times, preferably between 
7.5 and 11 times the nominal resonance period of 
the casing. The nominal resonance period of the 

20 casing is determined from the nominal thickness of 
the casing and the sound velocity In the casing, 
both parameters being known. 

For the window signals C(t) and P{t), Hamming 
windows can be used. Figure 4 shows such a 

25 Hamming window for C(t). centered on time origin 
to. As regards the processing window P(t). although 
Rgure 5 shows for it a symmetrical curve, it is 
preferred to use a Hamming window centered on 
the time origin to but truncated at time ti i.e. an 

30 asymmetrical window signal. 

Characterizing the resonance of the casing 

35 The resonance of the casing in response to the 

acoustic ' pulse depends on the cement-casing 
bond, which will be expressed as an acoustic im- 
pedance of the cement, and the casing thickness, 
which are the primary parameters of Interest In the 

40 casing inspection. These parameters can be deter- 
mined from a characterization of a resonance mode 
of the casing e.g. the fundamental mode. 

As noted above, the wide processing window 
P(t) is defined so as to contain the direct reflection 

45 from the casing and the early part of the reso- 
nance. The narrow window C(t) only contains the 
reflection from the casing and as such, the win- 
dowed signal CW(t) provides information on the 
system formed by the pulse generating assembly 

50 transducer and drive electronics - and the mud. In 
other words, it provides the response of the 
"Infinite-block" (i.e. a casing of Infinite thickness 
and extent). In the determination of the parameters 
of interest, the information contained In the signal 

55 portion CW will be used to normalize the informa- 
tion contained in signal S(t) within the processing 
window P(t) so as to automatically compensate for. 
variations (with temperature, time etc) in the char- 
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acteristics of the acoustic pulse and the attenuation 
of the mud. 

In order to find and characterize the resonance, 
the frequency spectra of the windowed signals PW 
and CW are used. Figure 6 shows in solid lines the 
amplitude spectrum PW(a») obtained from signal 
PW(t). using a discrete Fourier transform, and In 
dashed lines the amplitude spectrum CW(a)) ob- 
tained likewise from signal CW(t). The resonance of 
the casing appears as a marked dip in signal PW- 
(«). 

Now. the amplitude spectrum being highly de- 
pendent on transducer variations and attenuation 
by the borehole fluid, it is preferred to use the 
group delay spectrum. The group delay is defined 
as the derivative of radian phase vs. angular fre- 
quency and is computed from the phase spectrum. 
The group delay response is advantageous in that 
it tends to be smooth and linear except at resonan- 
ces and therefore, allows easy identification and 
characterization of the resonance. 

This is illustrated in Figure 8, which shows the 
theoretical amplitude (solid line) and group delay 
(dashed line) responses of a casing of 1 78 millime- 
ters (7 inches) diameter and 4 millimeters thick- 
ness, in the case of a solid cement bond (cement 
impedance of 8 MRayl). Figure 9 shows likewise 
the amplitude and group delay responses of a 
casing of 178 millimeters diameter and 16 mm 
thickness, with the same cement bond. 

Rgure 7 highlights the group delay spectrum 
around a resonance. The resonance can be char- 
acterized by the central frequency Fo of the reso- 
nance, corresponding to a minimum group delay 
rmin, the depth H of the dip corresponding to the 
resonance, and- the fractional bandwidth B/Fo. The 
bandwidth B is the width of the dip at a level At 
above the group delay minimum Tm,„ equal to a 
predetermined percentage, e.g. 40%, of the depth 
H of the dip : B = Fa - Fi , with the corresponding 
group delays ra and n such that t2 - T,„in = ti - 

rmm = At = 0.4H. 

The central frequency Fo is primarily depen- 
dent on the casing thickness while the depth and 
fractional bandwidth are primarily influenced by the 
acoustic impedances of the mud and the cement. 



Determining casing thickness and cement imped- 
ance — 

A normalization step of the group delay spec- 
trum is first performed In order to eliminate any 
slope, by linearizing the group delay spectrum of 
the normalization signal CW and subtracting it from 
the spectrum of the signal PW. 

Then the parameters of the resonance are de- 
termined, by searching a minimum in the normal- 



ized group delay spectrum within a pre-determined 
frequency range located around the nominal reso- 
nance frequency of the casing. This yields values 
for the parameters FO, H and B/FO. which will be 

s used as the '^measured parameters" in the follow- 
ing description. 

After the measured parameters of the reso- 
nance have been determined from the normalized 
group delay spectrum, the casing thickness and 

w cement impedance ,are determined by iteration us- 
ing a simple model, shown in Figure 10. This 
model is a planar model (i.e. the casing wail is 
modelled as a plane) defined by the cement im- 
pedance Zee and the casing thickness d. with mud 

75 impedance Zm and casing impedance Zp. Mud 
impedance Z^ is obtained by separate calibration 
and casing impedance Zp is a known data. Resis- 
tances R can be added to the model to perform 
corrections taking into account the curvature of the 

20 casing, the resistances being set to values defined 
in accordance with the diameter of the casing. 

The response of the model to a Dirac pulse e 
is determined in the frequency domain, with the 
initial value for the cement impedance Zee obtained 

25 from the measured depth H of resonance by using 
look-up tables, and the initial value of casing thick- 
ness derived from central frequency Fo. This yields 
an impulse response spectrum of the model. 

This spectrum is then multiplied by the nor- 

30 malizatlon spectrum CW(u) representative of the 
transducer-mud system to obtain what can be 
termed the transducer-mud spectrum. 

The processing window is then applied in the 
frequency domain by convolution with the 

35 transducer-mud spectrum, and the group delay 
spectrum of the convoluted signal is computed. 

This group delay spectrum is again normalized 
as described above by subtracting the linearized 
group delay spectrum of the normalization signal 

40 from the group delay spectrum of the convoluted 
signal. 

The normalized group delay spectrum is then 
used to determine "model" values for the param- 
eters of the resonance, namely the central fre- 

45 quency F„,od corresponding to the group delay 
minimum Tmid identified in this spectrum within the 
specified frequency range, and the fractional band- 
width B,„od taken at a level T,„od = ^-mid + At (At 
being as defined above). These model parameters 

50 are compared with the measured parameters. From 
the comparison, new estimates for the cement im- 
pedance and casing thickness are defined for a 
second iteration. The iteration is continued until the 
model group delay response is characterized by 

55 the same parameters as the measured resonance. 
Iteration may be performed by Newton-Raphson 
fitting technique, which typically requires only three 
iterations. 
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The planar model used in this method is quite 
simple and has the advantage of minimizing pro- 
cessing time. 

A more exact model, based on the cylindrical 
geometry of the casing, could be used instead of 
the planar model. Such a cylindrical model would 
include additional parameters (casing diameter, 
transducer dimensions, transducer offset from the 
borehole axis, sound velocity in mud). 



Alternative computation methods 

Other computation methods can be envisioned 
with the same basic definitions of time windows C- 
(t) and P(t) as specified above. An alternative meth- 
od consists of the following steps. 

The entire signal S(t) is converted into the 
frequency domain before any windowing. The 
spectrum Situ) thus obtained is normalized by the 
spectrum CW(a)) derived from the calibration signal 
CW(t). the normalized amplitude being the ratio of 
the amplitudes and the normalized phase the dif* 
ference between the phases. 

The processing window Is then applied to the 
normalized spectrum in the frequency domain, by 
convolution of the normalized spectrum with the 
spectrum of the processing window. 

The group delay spectrum is then calculated 
from the result of the convolution. This allows the 
resonance to be characterized as described above, 
thus yielding central frequency, depth and frac- 
tional bandwidth. 

From this point, it is possible to determine the 
cement impedance and casing thickness by direct 
computation, using pre-established look-up tables 
relating these parameters with the parameters 
which characterize the resonance. Another ap- 
proach is to use the model as defined above. The 
computation steps will be as explained above, ex- 
cept that the step of normalizing the group delay 
spectrum can be dispensed with. 

Claims 

1 . A method for determining at least one char- 
acteristic of a casing cemented in a borehole, such 
as cement bond and casing thickness, from a re- 
flected acoustic signal S(t) obtained by directing an 
acoustic pulse at a substantially normal incidence 
towards a radial sector of the wall of the casing, 
said pulse stimulating thickness resonance within 
the walls of the casing, comprising the steps of : 
defining a first time window corresponding to a first 
portion of signal S(t) including the initial reflection 
from the casing and subsequent acoustic returns 
due to resonance, defining a second time window 
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corresponding to a second portion of signal S{t) 
only including the initial reflection from the casing, 
and determining said casing characteristic from 
information related to resonance contained in said 
5 first time window while normalizing said information 
by information contained within said second time 
window. 

2. A method according to claim 1 , wherein said 
second time window is defined so as to include the 

10 amplitude maximunrt of signal S(t). 

3. A method according to claim 2, wherein the 
definition of the first and second time windows 
comprises the steps of determining the instant at 
which the amplitude maximum occurs in the signal. 

IS said instant being taken as a time origin to. and 
defining first and second time windows on the 
basis of this time origin. 

4. A method according to claim 3, wherein, for 
determining the time origin, an analytic signal Sa(t) 

20 of the response S(t) is determined: the amplitude 
maximum of the analytic signal is determined; and 
the time origin to is determined as being that in- 
stant at which the annplitude maximum occurs in 
the analytic signal Sa(t) 

25 5. A method according to claim 4, wherein the 

analytic signal Sa(t) is obtained by a Hilbert trans- 
form. 

6. A method according to claim 4, wherein the 
analytic signal Sa(t) is obtained by performing a 

30 Fourier transform, filtering to eliminate negative fre- 
quencies from the spectrum, and performing a 
further Fourier transform on the signal as obtained 
after filtering. 

7. A method according to claim 3. wherein the 
35 time origin is defined by performing the following 

steps: extracting the absolute values of the am- 
plitude maximum in the acoustic signal and of two 
extreme values situated on either side of said maxi- 
mum; determining a second degree polynomial fit- 
40 ting these three values; and defining the time origin 
as the instant at which the polynomial is at a 
maximum. 

8. A method according to claim 3, wherein the 
time intervals spanned by the time windows are 

45 defined by multiples of the nominal resonant period 
of the casing. 

9. A method according to claim 8. wherein the 
second time window is centered on the time origin 
to, and extends on either side of the time origin 

so over a period of time lying in the range from 1.5 
and 3 times the nominal resonance period of the 
casing. 

10. A method according to claim 8 or 9. 
wherein the first time window extends before the 

55 time origin to over 1.5 to 3 times the nominal 
resonance period of the casing and after the time 
origin over 6 to 12 times said nominal resonance 
period. 

7 
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11. A method according to claim 1, wherein 
said characteristic of the casing is determined by 
combining information derived from the frequency 
spectrum of the acoustic signal S(t) within the first 
time window and information derived from frequen- s 
cy spectrum of the acoustic signal S(t) within the 
second time window. 

12. A method according to claim 1 1 comprising 
the step of determining the frequency spectra of 

the first portion of signal S(t) and of the second to 
portion of signal S(t), respectively, normalizing the 
spectrum of the first portion by the spectrum of the 
second portion, and characterizing from said nor- 
malized spectrum the resonance of the casing, the 
parameters of the resonance being indicative of the is 
characteristics of the casing. 

13. A method according to claim 12. wherein 
the frequency spectra are the group delay fre- 
quency spectra. 

14. A method for acoustic inspection of a cas- 20 
ing cemented in a borehole, comprising the steps 

of obtaining a reflected acoustic signal S(t) by 
directing an coustic pulse at a substantially normal 
incidence towards a radial sector of the wail of the 
casing, said pulse stimulating thickness resonance * 25 
within the walls of the casing, defining a first time 
window corresponding to a first portion of signal S- 
(t) including the initial reflection from the casing 
artd subsequent acoustic returns due to resonance, 
defining a second time window corresponding to a 30 
second portion of signal S(t) only including the 
initial reflection from the casing, and determining at 
least one characteristic related to the casing from 
information related to resonance contained in said 
first time window while normalizing said information 35 
by information contained within said second time 
window. 

15. An apparatus for acoustic inspection of a 
casing cemented in a borehole, comprising a 
sonde adapted for displacement In the borehole 40 
which includes transducer means for obtaining a 
reflected acoustic signal S(t) by directing an acous- 
tic pulse at a substantially normal incidence to- 
wards a radial sector of the wall of the casing, said 
pulse stimulating thickness resonance within the 45 
walls of the casing; and signal processing means 

for performing the steps of defining a first time 
window corresponding to a first portion of signal S- 
(t) including the initial reflection from the casing 
and subsequent acoustic returns due to resonance. 50 
defining a second time window corresponding to a 
second portion of signal S(t) only including the 
initial reflection from the casing, and determining at 
least one characteristic related to the casing from 
information related to resonance contained in said 55 
first time window while normalizing said information 
by information contained within said second time 
window. 
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® Method for determining at least one (Characteris- 
tic of a casing cemented in a borehole, such as 
cement bond and casing thickness, from a reflected 
acoustic signal S(t) obtained by directing an acoustic 
pulse at a substantially normal incidence towards a 
radial sector of the wall of the casing, said pulse 
stimulating thickness resonance within the walls of 
the casing. 

The method according to the invention com- 
prises the steps of : defining a first time window 
corresponding to a first portion of signal S(t) includ- 
ing the initial reflection from the casing and subse- 
quent acoustic returns due to resonance, defining a 
second time window corresponding to a second por- 
tion of signal S(t) only including the initial reflection 
from the casing, and determining said casing char- 
acteristic from information related to resonance con- 
tained in said first time window while normalizing 
said information by information contained within said 
second time window. 
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